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Abstract We investigate two-loop higher-order binding corrections to 
the fine structure, which contribute to the spin-dependent part of the 
Lamb shift. Our calculation focuses on the so-called "two-loop self- 
energy" involving two virtual closed photon loops. For bound states, 
this correction has proven to be notoriously difficult to evaluate. The 
calculation of the binding corrections to the bound-state two-loop self- 
energy is simplified by a separate treatment of hard and soft virtual 
photons. The two photon-energy scales are matched at the end of the 
calculation. We explain the significance of the mathematical methods 
employed in the calculation in a more general context, and present re- 
sults for the fine-structure difference of the two-loop self-energy through 
the order of a^. 
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1 Introduction 



Ultra-precise measurements in atomic systems represent today of the most stringent available 
tests of fundamental quantum theories and a means for the determination of fundamental phys- 
ical constants with unprecedented accuracy The theoretical description of the bound states 
at a level of accuracy which matches the current experimental precision, which has reached 1.8 
parts in 10^^ and whose accuracy is to be improved in the near future demands a thorough 
understanding of the bound state including - among other effects - the relativistic, one-loop, 
two-loop and higher-order radiative, recoil, radiative-recoil, and nuclear-size corrections [^,^. 

We focus here on radiative corrections, which can be described - for atomic systems with low 
nuclear charge number - by a nonanalytic expansion in powers of the three parameters (i) a 
(the fine-structure constant), (ii) the product Za {Z is the nuclear charge number), and (iii) the 
logarithm ln[(Za)~^]. The expansion in powers of a, which is the perturbation theory parameter 
in quantum electrodynamics (QED), corresponds to the number of loops in the diagrams. The 
bound-state effects are taken into account by the expansions in the two latter parameters. 
Higher-order terms in the expansions in powers of Za and ln[(Za)~^] are referred to as the 
"binding corrections". One of the historically most problematic sets of Feynman diagrams in 
the treatment of the Lamb shift for atomic systems has been the radiative correction due to two 
closed virtual-photon loops shown in fig. |^. 




Figure 1: Feynman diagrams representing the two-photon 
electron self-energy. The double line denotes the bound elec- 
tron propagator. The arrow of time is from right to left. 

Let us recall at this point that even the evaluation of higher-order binding corrections to the 
one-loop self-energy, which a "priori should represent a less involved calculational challenge, has 
represented a problem for analytic evaluations for over three decades j^-^. The energy shifts 
of the bound states due to the radiative corrections are conveniently expressed by expansion 
coefficients corresponding to the powers of Za and ln[(Za)~^]; the naming convention is that 
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the power of Zol and the power of the logarithm are indicated as indices to the analytic coeffi- 
cients [see also eq. (0)] below. Because the expansion in both the one-loop and the two-loop case 
starts with the fourth power of Za, the non-vanishing coefficients carry indices A].i and -B^i for 
the one- and two-loop cases, respectively (with A; > 4 - see for a comprehensive review). 

Logarithmic corrections with / > 1 can sometimes be inferred separately in a much simplified 
approach, e.g. by considering infrared divergent contributions to electron form factors. By con- 
trast, the higher-order non- logarithmic coefficients represent a considerable calculational chal- 



lenge. Realistically, i.e. with the help of current computer algebra systems |1C,11|, one can hope 
to evaluate non-logarithmic coefficients of sixth order in Za. Complete results for the one-loop 
higher-order correction Ago for S and P states have only been available recently [P, 12], 13]. Calcu- 



lational difficulties have by now precluded a successful evaluation of the corresponding coefficient 



i?60 foi' the two-loop effect. Ground- work for the evaluation of -Bgo was laid in |14|. Here, we 
are concerned with the evaluation of the fine-structure differences of the logarithmic and non- 
logarithmic coefficients -Bgi (where L = 0,1,2), i.e. with the nP3/2~?T'Pi/2 difference of these 
coefficients. 

Using natural Gaussian units (^ = c = eo = 1), as it is customary for the current type of 
calculation, we write the two-photon self-energy in the Za-expansion for P states in terms of 
-B-coefficients as 



~3 



A^SE = - {ZaT — 



B. 



40 



+ {ZaY Bq2 \n\Zay^ + B^i In(Za)"^ + B, 



60 



+ 7^ 



where the remainder TZ is of order 0{ZaY . Relevant Feynman diagrams are shown in fig. ||. 

Here, m denotes the electron mass (we write eq. (|l]) in the non-recoil limit, i.e. for an infi- 
nite nuclear mass). The double logarithmic -B62-coefficient is spin-independent, so that we have 
Afs-B62 = 0. In this paper, we evaluate the fine-structure differences 

Afs^ei = ^6l("'P3/2) - ^6l("'Pl/2) , 

AfeSeo = S6o(nP3/2) - S6o(nPi/2) . (2) 

Throughout the paper, we will follow the convention that AfgX = X{nV^i2) — X(nPx/2) denotes 
the "fine-structure part" of a given quantity X. For Afs-Bgi and Afgi^gO) we provide complete 
results. It is perhaps worth noting that two-loop self-energy effects for bound states have repre- 
sented a considerable challenge for theoretical evaluations. Our investigation represents a contin- 
uation of previous work on the two-loop problem (see e.g. [|l4HT7||). It is probably a triviality to 
express that technical difficulties in the calculation and its description in the following sections 
of the paper cannot be avoided. 

For the description of the self-energy radiative effects - mediated by hard virtual photons -, we 
use the modified Dirac hamiltonian 

i/J,'^^ = a- [p-eFi(A) A] +/3m + eFi(A)(/) + F2(A)-^ {i^i ■ E - 13 a ■ B) , (3) 

which approximately describes an electron subject to an external scalar potential (f) = 0(r) and 
an external vector potential A = A(r). This modified hamiltonian is still local in coordinate 
space. The Dirac matrices in (P) are to be understood in the standard (Dirac) representation [|l^ 
(in the sequel, we will also use the non-covariant notation (3 = and a* = J^j'^)- 
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The argument A of the electron form factors Fi and F2 in eq. (^) is to be interpreted as a 
Laplacian operator acting on all quantities to the right (but not on the wave function of the 
bound electron in evaluating H^^\^p)). In momentum space, the action of the hamiltonian 
is described by the convolution H^^tp {p') = J d^p/(27r)^ H^\p' —p) 'ip{p). The form factors 

- in momentum space - assume arguments according to the replacement A —q^ = —{p' — p)'^- 
In eq. (^), radiative corrections are taken into account in the sense of an effective theory via 
the inclusion of the on-shell form factors Fi and F2. Although the bound electron is not an 
on-shell particle, the modified hamiltonian (^) can still approximately account for significant 
radiative systems with low nuclear charge number Z. Of course, the hamiltonian (^) cannot 
offer a complete description of the bound electron. Recoil effects cannot be described by a 
one-particle equation in principle, and vacuum-polarization effects are not contained in eq. (^). 
However, the effective description of self-energy radiative corrections mediated by hard virtual 
photons given by eq. (|3|) will turn out to be useful in the context of the current investigation. 

Both of the form factors Fi and F2 entering in eq. @ are infrared divergent, but this divergence 
is cut off in a natural way at the atomic binding energy scale (Za)'^ m. The fact that on-shell 
form factors can describe radiative corrections to the fine structure - mediated by high-energy 
virtual photons - has been demonstrated explicitly in |jl^ . The modified Dirac hamiltonian (^) 
and the associated modified Dirac equation have been introduced - in the one-loop approxi- 
mation - in ch. 7 of |jl^ [see for example eqs. (7-77) and (7-103) ibid.]. The low-energy part 
of the calculation is carried out using nonrelativistic approximations in the spirit of the sim- 
plified treatment introduced in the previous one- and two-loop calculations |,|l2|^0,||]. This 
approach was inspired, in part, by various attempts to formulate simplified low-energy ("nonrel- 
ativistic") approximations to quantum electrodynamics ("NRQED"), see e.g. [21,2^. Both the 
high-energy and the low-energy contributions are matched at the separation scale e whose role 
in the calculation is illustrated by the mathematical model example discussed in app. 

In a two-loop calculation, either of the two virtual photons may have a high or low energy as 
compared to the separation scale e. A priori, this necessitates jl^ a separation of the calculation 
into three different contributions: (i) both photon energies large, (ii) one photon with a large and 
one with a small energy, and (iii) both photons with small energies. For the particular problem 
at hand (the fine-structure differences of -Bgi and Bqq), we are in the fortunate position that 
effects caused by hard virtual photons (i) are described by the modified Dirac hamiltonian (^), 
whereas the low-energy part discussed in sec. § below comprises both remaining contributions 
(ii) and (iii). 

This paper is organized as follows: Two-loop form factors entering in eq. (^) are analyzed in 
sec. ^. The calculation is split into two parts: the high-energy part discussed in sec. |3| and the 
low-energy part, which is treated along ideas introduced in in sec. |^. Results and conclusions 
are left to sec. |5[ 



2 Two-loop Form Factors 

In order to analyze the modified Dirac hamiltonian (|3|) through two-loop order, we first have 
to investigate certain expansion coefficients of the electronic Fi and F2 form factors which 
are thoroughly discussed in the seminal papers ||2^,|2^. For the momentum transfer q"^ which 
is the argument of the two functions Fi = Fi{q'^) and F2 = F2{q^), we use the convention 
q^ = q^q^ = (q^)"^ — . The variable t in 24| is given as t = q^. When we evaluate radiative 
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corrections to the binding Coulomb field which is mediated by space-like virtual photons, we 
have = —q^ because for = 0. We use the conventions (see eq. (1.2) in [^ ): 

^iW = l + E© F.{t) = Y.il) Ft\t)- (4) 

n=l ^ 71=1 ^ ^ 

One and two-loop effects are denoted by upper indices 2 and 4, respectively. This notation is mo- 
tivated by the observation that two-loop effects are of forth order in the quantum electrodynamic 
interaction Lagrangian — e ^ ^^A^ %[) (in the Furry picture, which is used for the description of 
bound states, the Coulomb interaction is taken out of the interaction Lagrangian). 

There are two different points of view regarding the choice of diagrams to be included in the 
two-loop form factors, depending on whether the self-energy vacuum polarization diagram fig. ^ 
is included in the calculation or not. We will discuss both cases and give results with and without 
the diagram shown in fig. |2| taken into account. 




Figure 2: Combined self-energy vacuum-polarization diagram (denoted 
"1/" in the text). 



First, we discuss results obtained for Fi including the combined self-energy vacuum polarization 
diagram. In this case, the known results for the slope F[{<d) and for -F2(0), through two- loop 
order, read as follows. From eq. (1.11) of [23|, we have: 



a 



m^F[{{)) = - 

IT 



+ 1^ 

TT 



1 



-- In I — 

m 



A 



1 



2 r 



4819 49 3 

_ ^ ^(2) + 3 ((2) In 2 - - C(3) 
5184 72 ^^ ^ ^ 



where the forth-order coefficient has the numerical value 

,2W(4) 



According to eq. (1.7) in p 



m^F[^ '{<d) = 0.469 941487460. 
the value of -F2(0), through two-loop order, reads 



1 a 



^2(0) = --+ - 



2 TT 



a 



TT 



197 1 3 ■ 

— + -C(2)-3C(2) In2 + -C(3) 



(5) 



(6) 



(7) 
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where the two-loop coefficient has the numerical value 

Ff^(O) = -0.328478965 579. 

We now turn to the discussion of the slope F2*^'*^(0). In view of eq. (1-20) of |^] (see also ||2 
we have (up to two-loop order) 



(8) 



a 



a 



TT 



\n^B{t)4'\t)+4''\t) 



where the coefficients J- are by definition infrared safe and 



t 



(9) 



(10) 



Equations (^) and (|T0|) uniquely determine the infrared divergent contribution to -F2''^^(0). An 
analytic expressions for J^2^\t), t spacelike, has recently been obtained |^] in terms of harmonic 
poly logarithms |27,2^. As a byproduct, an analytic expression for the the slope J-2'^\o) was 
found. The result reads 



m2Ff)(0) 
m2j-f)(0) 



6 ym 



2160 20 ^^ ^ 10^^ ^ 40^^ ^ 



(11) 



A numerical result for J^2'^\o), complementing the above analytic expression, can easily be 
derived in combining eq. (1.20), eq. (1.30), and eq. (3.2) in |2^, as will be explained in the 
sequel. The dispersion relation (1.30) in [p3] reads. 



Am 1 t 

^ t — Am? ~'~ TT i — Am^ 



P 



dt' t' - Am^ 



t'-t t' 



ImF2(t') 



fl2) 



4m2 



where P denotes the Cauchy principal value. Equation (12) applies also if we single out the two- 
loop effect and differentiate at zero momentum transfer, and we obtain for the slope F2^\o) the 
relation 



m^Ff)(0) = lFf(0) + ^P / dt' 



where F2^^^(0) is given in eq. (pi). The second term on the right-hand side, denoted by T, can be 
evaluated using the result for ImF2^^(x) presented in eq. (3.2) in [23|; it reads 



l!I^ImFf'(O = lFf(0)+T. (13) 



4m^ 



T 



dx 



(1 - xf 



Im Fi^Vx) = -- ln( — ) +0.030 740 507833(1). 



1 



A 



^ x{l + x) ^ ' ' 6 

Here, the last error is due to numerical integration, and use is made of the natural variable 

_ 1 - y/1 - Am'-^/t 
^ ~ 1 + ^1 - A-r^Jt ' 



(14) 



(15) 
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In combining the result of eq. (|8|) with eqs. ( p^ ) and (H), the result m? J^2^\o) = 
—0.051379 233 561(1) is obtained which is in agreement with (|ll|). 

Now we will provide results for the form factors obtained excluding the self-energy vacuum- 
polarization graph V shown in fig. ^. These results refer to the pure two-photon self-energy dia- 
grams shown in fig. |^. The two-loop self-energy diagrams independently form a gauge-invariant 
set. They represent a historically problematic correction, and are the main subject of our in- 
vestigation. The combined self-energy vacuum-polarization diagram, according to eqs. (1.9) and 
(1.10) in |24] - taking into account the subtracted dispersion relation (1.30) of [|2^ - leads to 



the following corrections: 

^'^(0) = + ^ C(2) = 0.031588 972474, 

F^^'^'^iO) = ^ - 2C(2) = 0.015 687421859, 

f}^^'^{0) = — - -C(2) = 0.000497506 323. (16) 
216 8 

For the pure self-energy graphs, which we would like to denote by the symbol S, we therefore 
obtain the following results. 



m2Ff)'^(0) = - 11^ C(2)+3C(2)ln2-^C(3) = 0.438352 514986, (17) 



F2^^)'^(0) = -— + -C(2)-3C(2) In2 + -C(3) = -0.344166387438, (18) 
16 2 4 



-- In(-) -0.051876 739885 = -- In f- ) + X^^)'^(O) . (19) 
o \m J 6 \m' 



where the latter equality defines ^'•^^''^(0) in analogy with eqs. (^ and (|Tl]). 



3 High— Energy Part 

Based on the modified Dirac hamiltonian (^, corrections to the energy of the bound Dirac 
particle can be inferred. We will refer to the energy corrections attributable to the Fi and F2 
form factors as Ei and E2, respectively. For Ei, we have 



El 



Fi{-q^)-l 



(20) 



where the index fs refers to the fine-structure terms, i.e. to the result obtained by subtracting 
the value of the matrix element for a nP^j^ state from the value of the same matrix element 
evaluated on a nPi/2 state. A matrix element (^)fs of a given operator A is evaluated as 



V', 



nP 



3/2 



V'nP 



3/2 



nP 



1/2 



V'nPi 



/2 
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where denotes the hermitian conjugate of the Dirac wave function ip (not the Dirac adjoint 
ij) = -0+^'^). The Dirac wave functions ijj are expanded in powers of (Za) up to the order 
relevant for the current investigation. This expansion avoids potential problems associated with 
the logarithmic divergence of the Dirac wave function at the origin. 

For El, up to the order of (Za)^, we have 



(21) 



For P states, the nonrelativistic (Schrodinger) wave function - the leading term in the Za- 
expansion of the Dirac wave function ~ vanishes at r = 0, but the first relativistic correction 
gives a finite contribution, resulting in 



n 



1 



4n5 



(Zafm^ 



(22) 



This leads - again up to the order of (Za)^ - to the following result for Ei, 



El 



(Zaf 



Ff)(0) 



n 



m 



(23) 



Observe that the derivative of the Fi form factor has a physical dimension of in natural 

units, giving the correct physical dimension for Ei. The correction due to F2 in (^ reads, 



E2 = {F2{-q^)—ifE)i,. 



(24) 



A particle in an external binding Coulomb field feels an electric field E = \[Ze)q/q^ - in 
momentum space -01 E = —{Ze) r/(47rr'^) in coordinate space. Vacuum polarization corrections 
to E = — (Ze) r/(47rr^) lead to higher-order effects. The correction E2 splits up in a natural 
way into two contributions E2a and £'26 which are associated with ^2(0) and the slope -^2(0)' 
respectively. E2a reads 



E- 



2. = f^^Ff(0)(-.^ 

Zm \ r'^ 



The evaluation of the matrix element leads to 



. 7 • T 



+ 



fs 



23 



487 5 
360 4^ ~ 10n2 



{Zaf 



m 



(25) 



(26) 



For the purpose of the current investigation, the (Za)^-component of E2a is selected only: 



E- 



2a 



n{4), 



, /487 5 

= <°> (720 + 



23 \ 



m . 



(27) 



The matrix element E2b can be expressed as 

E2b 



2m 

A transformation into coordinate space leads to 

d 



(7-9) 



dx 



n 



1 



- x=0,fs 



{Zafm^ 



(28) 



(29) 
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As a function of the principal quantum number n, the result for E2b reads: 



E. 



2b 



{Zaf 



-2Ff)(0) 



n 



m 



(30) 



This result involves the infrared divergent slope of the F2 form factor [see eqs. ([Tl| ) and ([l9|)]. 
We are thus faced with the problem of matching the infrared divergence of the slope of the F2 
form factor, expressed in terms of the fictitious photon mass A, with the usual (energy matching 
parameter) e introduced originally in Q. This can be done in two ways: (i) by matching the 
infrared divergence of the rate of soft bremsstrahlung, calculated with a fictitious photon mass A, 
to a result of the same calculation, carried out with an explicit infrared cut-off e for the photon 
energy. This way of calculation is described on pp. 361-362 of It leads to the result 



In 



A_ 

Ye 



(31) 



The matching procedure (ii) consists in a comparison of the result of the application of the 
formalism considered above, and its application to the high-energy part of the ground state 
Lamb shift, which is in leading order given by the infrared divergence of the F\ form factor, and 
the result obtained by direct calculation of this high-energy part in a non-covariant formalism 
with an explicit energy cut-off e, as it has been carried out in Q. This second matching procedure 
leads to the following result - in agreement with (|3l|) -, 



In- 



So, we are led to the replacement 



In- 



, m 11 

In \ 

2e 24 



m 5 

In h - 

2e 6 



(32) 



(33) 



A comparison with the results in eqs. (pT|), (^9|), and (| 
for the fine-structure difference is given by a term 



reveals that the logarithmic divergence 



n 



1 , m 

— IT- In — 

3n2 2e 



so that we may anticipate at this stage the result for AfsSgi, 

^2 1 



n 



fs-D61 



3n2 



Based on (|3^) and (| 

F2h - 



we can express E^h in terms of e and .7^*-^''(0), 



{Zaf 



1 

3 



In^- A + 2^f) 
2e V18 2 



(O)m^ 



(34) 



(35) 



n 



m . 



(36) 



There is a third correction due to the effect of two one-loop corrections on the electron vertices. 
Because we are only interested in the fine structure, we isolate the terms which are proportional 
to the spin-orbit coupling, and we obtain 

E3 = (^mm [-^^J mm f^^^^ , (37) 
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where 



Za 



cr ■ L , 



(38) 



and 1/{E — H)' is the nonrelativistic, spin-independent reduced Schrodinger-Coulomb Green 
function |^,|3^. The only spin-dependence in (p7|) occurs in the coupling a ■ L, and it can be 
taken into account by an overall factor, 



(39) 



We are therefore led to consider the "spin-independent version" of the matrix element which 
occurs in eq. ( |37[ ) and obtain the following result. 



Za 



1 



Za 



nP 



227 
8640 



1 1 
96n ^ 80n2 



[Za)^ m 



(40) 



The spin-dependence can be easily restored by considering eq. (|39D. The index "nP" in 
eq. (^0|) means that the matrix element is evaluated with the nonrelativistic, spin- independent 
(Schrodinger) wave function. Alternatively, on may evaluate with either the nPi/2 or the nP3y2 
Dirac wave function and expand up to the leading order in (Za). 

The evaluation of (^0|) can proceed e.g. by solving the differential equation which defines the 
correction to the wave function induced by Hfs, and subsequent direct evaluation of the resulting 
standard integrals using computer algebra [llO, 11|. The final result for E-^ reads. 



E; 



(Za) 



6 r 



227 1 
2880 ^ m 



80n2 



m . 



(41) 



This concludes the discussion of the high-energy part. The final result for the high-energy part 
is 

E^ = E, + E2a + E2b + Es , (42) 
where E'l, E2a, E2b, E^ are given in eqs. (p3), (p7|), (p6|), (41), respectively. 



4 Low— Energy Part 

The low-energy part consists essentially of two contributions. Both effects, denoted here -^4 
and E^, can be obtained by a suitable variation of the low-energy part of the one-loop self 
energy, by considering the spin-dependent effects introduced by a further one-loop electron 
anomalous magnetic moment interaction. The first of the two terms, E/^, is caused by spin- 
dependent higher-order effects in the one-loop self-energy, which receive additional corrections 
due to the anomalous magnetic moment of the electron. The second term, E^, is due to an 
anomalous magnetic moment correction to the electron transition current, which can also be 
seen as a correction to the radiation field of the electron due to its anomalous magnetic moment. 

The leading-order low-energy part (see ^) reads 

1 

H-{E-oj) 

In order to isolate the fine-structure effects, we should now consider corrections to the wave 
function, to the current, to the hamiltonian and to the energy of the bound state due to the 
spin-dependent relativistic (spin-orbit) hamiltonian 

^ / ^ e „ a(Za) —iy-r , , 

n = F2{0)—ij-E = ^ 44 

2m Airm r-^ 



2a r , 

El = / dujuj 

Sirm Jo 



(43) 
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The above hamiltonian TC is the last term in the modified Dirac hamiltonian [right-hand side 
of eq. @], approximated for a particle bound in a Coulomb field with the F2 form factor 



evaluated at zero momentum. The electric field E in (44) corresponds to the binding Coulomb 
interaction. The hamiltonian describes the modification of the spin-orbit interaction due to 
the anomalous magnetic moment of the electron. 

multiplied by 
The resulting 



The nonrelativistic limit of Tl is the spin-orbit coupling Hf^ given in eq. ( |38|) 
a factor 2i<2^^(0) = q/tt (the additional factor 2 finds an explanation in 
hamiltonian 



H. 



a 



a 



Za 



cS 



^fe = - -T-Y^ ^L. (45) 

takes into account magnetic vertex corrections in the framework of an effective theory. Denoting 
the variation of the expression ( ^3|) mediated by Hcs with the symbol 6es ~ in the spirit of the 
notation introduced in [141 -, we obtain the contribution 



4a 



5ofr 



2q 
Svrm 



dw w ( (/) 



1 



H-{E-cj) 



P 



(46) 



Following the notation introduced in [fL2| , [l^ , the contribution E4a is the sum of the fine-structure 
effects created by the wave- function-correction Fs^, the first relativistic correction to the energy 
FsE, and the correction due to the relativistic hamiltonian Fsh, each multiplied by a factor a/vr. 
That is to say, the final result for E4 is 



(47) 



vr 



n 



There is a further correction to the nonrelativistic effective coupling to the radiation field due 
to the "anomalous spin-orbit hamiltonian" (^). The correction, in the nonrelativistic limit, can 
be derived by considering a Foldy-Wouthuysen transformation which by definition diagonalizes 
the hamiltonian (|4^ in spinor space and also leads to a transformation of the relativistic current 
operator a* according to 



a 



Ua'U- 



U 



exp 



2m 



Here, (5 and a* are standard Dirac matrices [19[, i is a spatial index, and H is given in 
calculation is carried out along ideas introduced in [02] and leads to the result 



5j 



a Za 



46 



vr 



a X r . 



(48) 



The 



(49) 



as a relativistic correction to the electron current which is simply a* in the relativistic formalism 
and p'/m in the leading nonrelativistic approximation. Again, following the notation introduced 
in [12,13|, the resulting additional contribution is 

E,b=(-y {Za)^^Ai,Fsy. (50) 
The sum of (^) and (^) is just the (Za)®-component of the fine-structure difference of the 



one-loop self energy from [12, 13[, multiplied by an additional factor a/vr. It can also be written 

as 



E4 — Eia + E41, 



(Za) m 



n 



1 



3^2 



In 



2e 



+ 



n 



1 



[ZaYm n 



(51) 
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where Afs^4(n) could be interpreted as a relativistic generalization of a Bethe logarithm, which 
is n-dependent. However, a significant numerical fraction of the n-dependence can be eliminated 
if the factor (n^ — l)/n^ is taken out of the final result. The evaluation of Afsii{n) has recently 
been performed in |32] with improved numerical methods (see e.g. [Q), and the following results 
have been obtained: 

AM2) = 0.512 559 768(1), 
AfshiS) = 0.511 978 815(1) , 
Afj4{4) = 0.516 095 539(1), 

Afe^4(5) = 0.519 976 941(1) , (52) 
where the uncertainty is due to numerical integration. 

There is, as stated above, a further correction due to the explicit modification of the transition 
current due to the anomalous magnetic moment; it can be obtained through the replacement 



a' a* + F2(0) 



2m 



(53) 



and must be considered in addition to the correction (p8[). A careful consideration of the non- 
relativistic limit of this correction to the current, including retardation effects, leads to the 
result 

r • a Za , ^, 

oj^ = 5- a X r . (54) 

•^^ 27r 2mr3 ^ ^ 

Consequently, we find that the correction is effectively -^2(0) times the retardation corrections 
to the transition current Fsy found in [12,1^. We obtain 



E5 



a 



4 m AfsFsy 



,7r/ ' ' 2 
In analogy with E4, this correction can favorably be rewritten as 



n 



1 



Afe4(n) 



(55) 



(56) 



On the basis of ||T2, T^, 3^], we obtain 



Afs4(2) 
Afs4(3) 
Afe4(4) 
Afs4(5) 



-0.173 344 868(1) 
-0.164 776 514(1) 
-0.162 263 216(1) 
-0.161 165 602(1) 



The final result for the low-energy part is 

-E'l = -E4 -I- £'5 , 



(57) 



(58) 



with £'4 and £5 being given in eqs. ( [5l| ) and (5(:), respectively. 

We can now understand why it was possible to join the two contributions with "mixed" and "low- 
and-low" energy virtual photons (ii) and (iii), which were discussed in sec. |l], into a joint "low- 
energy part". The reason is simple: The effective hamiltonian (^) has no infrared divergence, 
because it involves the low-energy limit of the magnetic form factor £2, which is infrared safe 



12 



in one- loop order according to eq. (P). Because the main contribution to the quantity ^2(0) is 
caused by hard virtual photons, it is also justified to say that the contribution of "low-and-low" 
energy virtual photons vanishes at the order of interest for the current calculation (fine-structure 
difference). In higher- loop order, the further infrared divergence acquired by F2 would lead to 
an infrared divergence in the effective hamiltonian constructed in analogy with eq. ([45|); this 
infrared divergence would have to be attributed to a "mixed" contribution (one photon of high 
energy, and one low-energy photon). 



5 Results and Conclusions 

We have obtained analytic results for higher-order correction to the two-loop self-energy of P 
states in hydrogen-like systems. In our calculation, we have analyzed the electron form factors 
through two-loop order in sec. and we have split the calculation into a high-energy part with 
two hard virtual photons discussed in sec. ^, and a low-energy part with at least one soft virtual 
photon analyzed in sec. |^. The final result for the contribution to the fine-structure energy 
difference is obtained by adding the high-energy contributions Ei - given in eqs. (|23|), (|27|), 
dH), (|4l|), and the low -energy effects £'4 and E^ from eqs. ( [5l] ) and (|56|). The dependence on e 
cancels out in the final result which is the sum of the high-energy part En given in eq. (^) and 
the low-energy part E'l defined in eq. (|58[). This is also evident when considering explicitly the 
eqs. ( ^6|) and (|5l|). The final results for the analytic coefficients of order a'^{Za)^ read 

- 1 

AfsS6i = -^^- (59) 



[see also eq. (pSj)] and 



f 227 1 3 \ (4)s. /487 5 23 \ 
- I288O + m " 80^? j + ^2 ' [n6^M~^^) 



n2-l 



(60) 
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where explicit numerical results for F[^^^'^ (0), F2^^'^ (0) and ^2^^'*''^(0) can be found in eqs. (p^), 
([1^ and (|l9|), respectively. This result refers to the pure self-energy diagrams in fig. I. The 
result reads numerically for the principal quantum numbers n = 2-5, 

Afs56o(2) = -0.361 196 470(1) , (61) 

Afs56o(3) = -0.411 156 068(1) , (62) 

Afs56o(4) = -0.419 926 624(1), (63) 

Afs56o(5) = -0.419 832 876(1) . (64) 

If it is desired to add in the combined self-energy vacuum-polarization diagram from fig. ^, then 
the form-factor results from eqs. (^), (^) and (^) instead of the pure self-energy results given 



in eqs. ([17D, dig ) and (19) have to be used in evaluating (pOj). When including the combined 
self-energy vacuum-polarization diagram from fig. ^, there is no further low-energy contribution, 
so that the alternative set of numerical values for the form factors from eqs. (^, (^) and ( [Tl| ) 
fully takes into account the additional effect of the diagram in fig. |l] on the fine-structure in the 
order of (Za)^. 
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It is perhaps worth mentioning that for the one-loop self-energy, analytic coefficients are known 
only up to the order of a{Za)^ ||l^, but the remaining uncertainty is removed by recent non- 
perturbative numerical calculations 32,^. For the two-loop effect, the (Za)-expansion con- 



verges more rapidly than for the one-loop effect in absolute frequency units because of the 
additional radiative factor a/ir which decreases the overall size of the effect. 

It is hoped that the analytic calculations for low nuclear charge number Z will be supplemented 
in the future by an accurate numerical treatment of the two-loop self-energy problem (see also 



related recent work in the high-Z region, refs. [35-3^). This presupposes that the considerable 



numerical problems in the domain of small nuclear charge could be solved by adequate numerical 
methods, and that the further problem of the increased computational demand of the two-loop 
effect in comparison to the one-loop problem 34 1 can be tackled - possibly by massively 



parallel computer architectures. Note, however, that the most accurate theoretical predictions 
could only be reached in combining numerical and analytic results. The reason is the following: 
All numerical calculations are performed in the non-recoil limit which is the limit of infinite 
nuclear mass. This is not quite sufficient for an accurate theoretical treatment because the 
self-energy effect for a bound-state depends genuinely on the ratio of the orbiting particle to 
the nuclear mass - an effect beyond the recoil correction. For example, the argument of the 
logarithms in (||) should be replaced according to ln[(Za)~^] ln[a where a = m/irij: 

and rrir is the reduced mass The possibility to include these tiny, but important effects 
depends crucially on a reliable knowledge of the analytic coefficients in combination with an 
accurate numerical treatment of the problem. 

The analytic results can be used to obtain improved theoretical predictions for the hydrogenic 



fine structure as compared to the previous order-a ''-calculations |12,1S], because they remove 
the principal theoretical uncertainty in the order of due to the problematic two-loop self- 
energy which is represented diagrammatically in fig. |^. A compilation of the other corrections 
relevant at the order of a^, including but not limited to the vacuum polarization effects, whose 
evaluation is rather straightforward, will be presented elsewhere. Our calculation illustrates the 
usefulness of the simplified effective treatments of two- loop effects in the analytic approach based 
on the modified Dirac hamiltonian (^) and the "e method" (see [§,1^, 1^, 20 1 and app. ^). This 
aspect highlights, as we believe, the need for systematic, simplified treatments of higher-order 
radiative corrections in bound systems. 

In this paper, we primarily address spin-dependent effects in one-electron (hydrogenlike) sys- 
tems. However, the same effects also contribute to the fine-structure splitting in two-electron 
(heliumlike) systems. There is currently remarkable interest in improved measurements of the 
fine-structure splitting in helium and heliumlike atomic systems with low nuclear charge 
The effects addressed in this paper contribute to the fine-structure splitting in helium on the 
level of 100 Hz, which is not much smaller than the current experimental accuracy of about 
1 kHz, and allows for an estimate of uncalculated yet higher-order contributions. 

The results of this paper are a step in a systematic study of the higher-order binding corrections 
to the two-loop Lamb shift of S and P states. The scheme of calculation permits not the only a 
simplified treatment of the problem via a separation into appropriate energy regions for the two 
virtual photons, but also a clear identification of the spin-independent and the spin-dependent 
contributions to the self-energy. The results will therefore be directly applicable to the total 
coefficients for P states once the spin-independent parts (the two-loop Bethe logarithms) are 
calculated. These are currently under study. 
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A The "e Method" 

We discuss here, by way of example, the e method employed in the analytic calculation of 



self-energy effects in bound systems. This method is very suitable [12, 13| for the separation 
of the two different energy scales for virtual photons: the nonrelativistic domain, in which the 
virtual photon assumes values of the order of the atomic binding energy, and the relativistic 
domain, in which the virtual photon assumes values of the order of the electron rest mass. 
Different approximation schemes and different asymptotic expansions are adequate for the two 
different domains. Without these approximations and expansions, the analytic evaluation of 
either the high- or the low-energy part would not be feasible. At the same time, the model 
example discussed in this appendix is meant to illustrate the usefulness of the "e method" in a 
more general context. 

We will consider here a model problem with only one "virtual photon". The separation into 
high-and low-energy photons necessitates the temporary introduction of a parameter e; the 
dependence on e cancels when the high- and the low-energy parts are added together. We have, 

nonrelativistic domain ^ e ^ electron rest mass (65) 
(Za)^me ^ e <C m-e , (66) 

where a is the fine structure constant, and Z is the nuclear charge. The high-energy part is 
associated with photon energies a; > e, and the low-energy part is associated with photon 
energies a; < e. 

In order to illustrate the procedure, we discuss a simple, one-dimensional example: the evaluation 
of 

■'■^ ~" dcj. (67) 

/o 

where the integration variable w might be interpreted as the "energy" of a "virtual photon". 
The integral / can be expressed in terms of special functions, 

/(/3)=;3E(-J,)=3|,F,(-l,i;l;-J,), (68) 

where E is the complete elliptic integral of the second kind, and 2-^1 denotes a hypergeometric 
function. An alternative integral representation reads /(/3) = /q^''^ \J ^"^ + sin^(u;) dw. 

The purpose of the calculation is to derive a semi-analytic expansion of I{0) in powers of (3 
and ln/3. The fine structure constant a takes the role of the expansion parameter [3 in actual 
self-energy calculations. We discuss first the "high-energy part" of the calculation. It is given by 
the expression 

W) = (69) 
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For uj > e, we may expand 



/32 /54 



^2 + ^2 = ^ + ^ + ^ + 0(/3^, (70) 

but this expansion is not applicable in higher orders to the domain < a; < e because of the 
appearance of inverse powers of uj (analogous to an "infrared divergence" in QED). 

The separation parameter e acts an infrared regulator. After expanding in /3 [see Eq. (|70|)], the 
resulting integrals in each order of (3 can be evaluated analytically. Subsequently, we expand 
every term in the /3-expansion in powers of e up to the order e'^, i.e. we keep only the divergent 
and constant terms in e. The result is 



/H(At) = l+f!'|ilii(|')+0(e)} 



+0{f3^) . (71) 

Here, the "C'-symbol identifies a contribution for which 0{x)/x const, as x ^ 0, whereas 
the "o"-symbol identifies the weaker requirement o(x) ^ as x — > 0; this is consistent with the 
standard notation (see e.g. |]38|| ). 

The contribution /h(/3) corresponds to the "high-energy part" in analytic self-energy calcula- 
tions, where the propagator of the bound electron is explicitly expanded in powers of the fine 
structure constant a. Now we turn to the "low-energy part". The expression for the low-energy 
part (0 < < e) reads 

HP) = Pxh^^du;. (72) 



The expansion (70) is not applicable in this energy domain; we therefore have to keep the numer- 



ator of the integrand y^uj^ + fS^ in unexpanded form. However, we can expand the denominator 
Vl — of the integrand in powers of uj; because < u; < e (with e small), this expansion in uj 
is in fact an expansion in (3 - although the situation is somewhat problematic in the sense that 
every term in the w-expansion gives rise to terms of arbitrarily high order in the /3-expansion 
[see also Eq. (^) below]. 

The term y^uj^ + (P is analogous to the Schrodinger-Coulomb propagator in the self-energy 
calculation which has to be kept in unexpanded form, whereas the expansion 

7T^-l + T + r' + °(-') (73) 
corresponds to the expansion into the (Za)-expansion in the low-energy part. 

Every term in the expansion ( [7^ ) gives rise to arbitrarily high-order corrections in (3, but it 
starts with the power uj"' — > For example, we have for the leading term of order o;^ = 1 

from Eq. (^, 

^^da; = jiln g e) + ^ + 0(e)} 

{ 16l^ + ^^'^] + {"61? + ^^'^] + ■ 
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Note that the terms generated in the orders and (3^ are needed to cancel divergent contribu- 
tions in respective orders of /3 from the high-energy part given in Eq. (70). The term of order 
o;^ from (72) results in 



^ Jo ^ 



1 



16 



/3 



1 



64 



+/3^{-^ + 0(6)}+0(/?«) 



(75) 



Altogether, we obtain for the low-energy part, 



Il(/3,6)=/32 jiln 



1 



1 

16^ 
1 



16 



In 



e + 



64 



\_ J_ /2 

64 £4 64 £2 + 128 ^\^^, 



+0(/3« ln/3) 



(76) 



When the high-energy part (|7lD and the low-energy part (|76|) are added, the dependence on e 
cancels, and we have 



m 



Inif3,e) + lUf3,e) 



2 V/3 



1 



+ 



+r -- In - + 777 



16 



13 



128 



In 



64 



-I 

128 J 



0(/3^ ln/3) 



(77) 



In order to illustrate the analogy with the self-energy calculation presented here, we would like 
to point out that the dependence on e cancels out in the final result which is the sum of the 
high-energy part -Eh given in eq. (42) and the low-energy part E'l in eq. (p8[). 
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